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We present a case that current observations may already indicate new gravitational physics on 
cosmological scales. The excess of power seen in the Lyman-a forest and small-scale CMB exper- 
iments, the anomalously large bulk flows seen both in peculiar velocity surveys and in kinetic SZ, 
and the higher ISW cross-correlation all indicate that structure may be more evolved than expected 
from ACDM. We argue that these observations find a natural explanation in models with infinite- 
volume (or, at least, cosmological-size) extra dimensions, where the graviton is a resonance with 
a tiny width. The longitudinal mode of the graviton mediates an extra scalar force which speeds 
up structure formation at late times, thereby accounting for the above anomalies. The required 
graviton Compton wavelength is relatively small compared to the present Hubble radius, of order 
300-600 Mpc. Moreover, with certain assumptions about the behavior of the longitudinal mode on 
super-Hubble scales, our modified gravity framework can also alleviate the tension with the low 
quadrupole and the peculiar vanishing of the CMB correlation function on large angular scales, 
seen both in CORE and WMAP. This relies on a novel mechanism that cancels a late-time ISW 
contribution against the primordial Sachs- Wolfe amplitude. 



I. INTRODUCTION 

In the ten years since the discovery of cosmic accelera- 
tion by the supcrnovae teams [1], the A-cold dark matter 
(ACDM) model has emerged as the standard paradigm 
for cosmology. With a single parameter (A), the model 
predicts expansion and growth histories consistent with 
observations. As near-future experiments will soon probe 
the large scale structure with unprecedented accuracy, it 
is worth asking at this juncture whether ten years from 
now the ACDM model will be hailed as the ultimate cos- 
mological theory, or whether it will have by then joined 
its Einstein de Sitter predecessor in the ranks of defunct 
cosmologies. A host of recent observations may have al- 
ready revealed cracks in ACDM's armor: 

1. The cross-correlation between galaxy surveys and 
the cosmic microwave background (CMB), result- 
ing from the Integrated Sachs- W'olfc (ISW^) effect, 
has recently been measured at the 4(7 level using the 
Wilkinson Microwave Anisotropy Probe (WMAP) 
data and a combination of large-scale structure sur- 
veys [2, 3]. This provides non-trivial and indepen- 
dent evidence for cosmic acceleration. Remarkably, 
the combined amplitude is larger than the ACDM 
prediction by > 2a's: 2.23 ± 0.60 [2], with 1 being 
the ACDM expectation. 

2. Independent recent analyses suggest anomalously 
large bulk flows on large scales. Using a com- 
pilation of peculiar velocity surveys, Watkins et 
al. [4] obtained a bulk flow of 407 ± 81 km s"^ 
on 50/i~^ Mpc scales, in contrast with the ACDM 
expectation of ~ 190 km s~^ for the r.m.s. value of 
the bulk flow, which are inconsistent at the 2(t level. 
Meanwhile, using the kinetic Sunyaev-Zeldovich ef- 
fect to estimate peculiar velocities of galaxy clus- 
ters, Kashlinsky et al. [5] found a coherent bulk 
motion of 600-1000 km s'^ out to > 300/1"! j^p^ 



3. High resolution observations of the CMB at i ^ 
3000 by the Cosmic Background Imager (CBI) in- 
dicate an excess power compared to the theoret- 
ical primary CMB power spectrum [6]. This ex- 
cess could be explained by the thermal Sunyaev- 
Zel'dovich effect from galaxy clusters at z ~ I. 
However, the amplitude of matter power spectrum, 
necessary to explain CBI excess is about 35% (at 
~ 2(7 level) larger than the current best-fit concor- 
dance ACDM model [7]. 

4. The Lyman-a forest in the spectra of high red- 
shift quasars provides us with the highest resolution 
precision measurement of the cosmological matter 
power spectrum. However, Lyman-a forest con- 
straints on the linear matter power spectrum, based 
on the SDSS quasar spectra at z ^ 3 [8, 9] is also 

35% (at ~ 2.5cr level) larger than the ACDM 
prediction (Fig. 1 in [10]). 

5. As first observed by the Cosmic Background 
Explorer (COBE) [11] and later conflrmed by 
WMAP [12], the CMB temperature anisotropy 
shows a lack of correlation on large (> 60 °) angular 
separations. While related to the low quadrupole, 
the real-space anomaly is more robust statistically 

less than 0.03% of random realizations exhibit 
smaller power on those scales [12, 13, 14]. The sta- 
tistical significance of this anomaly has only grown 
with subsequent WMAP data releases [15]. 

Of course, these anomalies may eventually go away, as 
systematic effects are better understood. After all, these 
represent for the most part ~ 2a discrepancies. But it 
is striking that the first four observations on this list of 
misfits all point towards the same physics: that structure 
is more evolved on large scales than predicted by ACDM. 
This suggests that gravity may be stronger at late times 
and on large scales. 



In this paper wc show that this is precisely what hap- 
pens in infrared- modified gravity theories, inspired by 
brane-world constructions with infinite-volume extra di- 
mensions. Because the extra dimensions are infinite in 
extent (at least cosmologically large), the 4d graviton is 
no longer massless but is instead broadened into a res- 
onance — a continuum of massive states — with a tiny 
width r~^. At first sight this may seem at odds with our 
earlier conclusion, since a massive/resonance graviton 
implies weaker gravity. While this is certainly true on the 
largest distances, on intermediate (albeit cosmologically- 
relevant) scales the longitudinal (or helicity-0) mode me- 
diates an extra scalax force which enhances gravitational 
attraction by order unity. 

This extra force is of course at the origin of the well- 
known van Dam-Veltman-Zakharov (vDVZ) discontinu- 
ity [16]. As conjectured by Vainshtein [17], however, non- 
linear interactions can suppress the effects of the longi- 
tudinal mode near astrophysical sources. This screening 
is also at play cosmologic;ally, suppressing the extra force 
at early times and on small scales. Thus gravity becomes 
stronger only at late times and on sufficiently large scales. 
This suppression at high density or curvature is qualita- 
tively similar to the chameleon mechanism [18, 19, 20], 
except that the Vainshtein effect arises from derivative 
interactions as opposed to a scalar potential. 

The theories of interest are higher-dimensional gen- 
eralizations of the Dvali-Gabadadze-Porrati (DGP) 
model [21], in which our visible universe is confined to 
a 3-brane. The confrontation of its cosmological predic- 
tions with observations has been the subject of much lit- 
erature since its advent [22, 23, 24, 25, 26, 27, 28, 29, 30]. 
Instead of having one extra dimension, here the bulk con- 
sists of 6 or more space-time dimensions. 

Extending the DGP scenario to higher dimensions has 
proven difficult historically. Early studies demonstrated 
that naive extensions are plagued by ghost-like instabili- 
ties, even around flat space [31, 32]. Recently, however, it 
was shown that such instabilities are absent if our 3-brane 
lies within a succession of higher-dimensional branes, 
each with their own induced gravity term, and embed- 
ded in one another in a flat bulk space-time [33, 34]. 
We refer to this framework as Cascading Gravity [35]. 
In the simplest codimension-2 case, for instance, our 3- 
brane is embedded in a 4-brane within a 6-dimensional 
bulk. A similar cascading behavior of the gravitational 
force law was also obtained in a different codimension-2 
framework [36]. 

Due to the higher-dimensional nature of these con- 
structions, extracting cosmological predictions presents a 
daunting technical challenge. Even deriving the modifled 
Friedmann equation is highly non-trivial. Here we cir- 
cumvent these difficulties by developing a phenomenolog- 
ical description for the background cosmology, while us- 
ing the parametrized post-Friedmann (PPF) [37] frame- 
work to encode modiflcations to Einstein gravity into 
the evolution of perturbations. This approach has been 
shown to accurately reproduce the cosmological behavior 



of the standard DGP model. Inspired by these results, we 
can infer the parametric dependence of the required PPF 
input functions and parameters for Cascading Gravity 
models. This allows us to calculate various cosmological 
observables, all within a 3-|-l-dimensional framework. 

One of our key observations is that in 6 or more 
dimensions the background evolution becomes nearly 
indistinguishable from ACDM. This allows us to con- 
sider relatively small graviton Compton wavelengths, e.g. 
Vc ~ 300 — 600 Mpc, without violating constraints on the 
expansion history. Compared to the usual assumption 
Tc ^ Hq^ in DGP, this expands considerably the dura- 
tion and range of scales over which enhanced gravity is 
effective. 

As mentioned earlier, the culprit for most of the inter- 
esting phenomenology in these models is the longitudinal 
or helicity-0 mode of the massive graviton. This degree 
of freedom is ultimately responsible for addressing all of 
the aforementioned anomalies: 

• The longitudinal degree of freedom modifies the 
time-evolution of the lensing potential, which re- 
sults in a stronger ISW cross-correlation compared 
to ACDM. 

• In the non-relativistic limit, the helicity-0 mode can 
be thought of as mediating an extra scalar force, 

which enhances gravitational attraction on inter- 
mediate scales. This speed-up in structure forma- 
tion at late times results in larger bulk flows, more 
clusters at 2; ^ 1, and enhanced power on Lyman-a 
scales. 

• As a consequence of this mode, our modified per- 
turbation equations exhibit an effective anisotropic 

stress component. With certain assumptions about 
the form of this component on super-Hubble scales, 
we find that the late-time ISW effect can destruc- 
tively interfere with the primordial Sachs-Wolfe 
amplitude, resulting in a power deficit in the CMB 
on large scales. 

Looking forward, our modified gravity models make 
several key predictions that will be tested by near-future 
experiments. For instance, because photons do not cou- 
ple directly to the longitudinal mode, the weak lensing 
potential is considerably less affected than the Newtonian 
potential. Thus a distinguishing prediction is a discrep- 
ancy in the value of erg estimated from the matter power 
spectrum and that inferred from either CMB or weak 
lensing observations. For Tc ^ 300 — 600 Mpc, the differ- 
ence is a 20-30% effect, depending on the redshift of the 
observation. 

The paper is organized as follows. In Sec. II, we de- 
scribe the building blocks of higher-dimensional gener- 
alizations to DGP, following [38, 39]. The result is a 
two-paramater family of massive or resonance gravity 

theories, with a specifying the form factor, and the 
Compton wavelength for the graviton. Motivated by the 
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DGP Fricdmann equation, wc propose in Sec. Ill a gen- 
eralized Friedmann equation for higher-dimensional the- 
ories. We then turn to cosmological perturbations in 
Sec. IV, and motivate our fiducial input parameters for 
PPF. Using these, we derive the resulting CMB power 
spectrum and angular correlation function in Sec. V. In 
Sec. VI, we study the predictions for structure formation, 
including the weak lensing power spectrum (Sec. VIA), 
the Integrated Sachs- Wolfe cross-correlation (Sec. VI B), 
large scale bulk fiows (Sec. VIC), the Lyman-a forest 
(Sec. VI D), and secondary CMB anisotropies from ther- 
mal Sunayev-Zeldovich (SZ) efi'ect (Sec. VIE). We con- 
clude with a summary of our results in Sec. VII. 



II. PHENOMENOLOGY OF 
HIGHER-DIMENSIONAL GRAVITY 



We are interested in generalizing the phenomenology 

of the DGP model [21], arguably the simplest example of 
modified gravity in the infrared. In this scenario, our visi- 
ble universe is confined to a 3-brane in a 4+1-dimensional 
bulk. Although the bulk has infinite volume. 4Z) gravity 
is recovered at short distances, due to an induced gravity 
term on the brane. The gravitational force law therefore 
scales as the usual at short distances, but weakens 
to l/r"^ at large distances, with the cross-over scale Vc set 
by the bulk and brane Planck masses: Vc = M^/M^. 

In momentum space, this corresponds to the graviton 
propagator l/(fc^ -|- r~^k). Following [38, 39], we will 
consider the following power-law generalizations: 



p^^-2(l-a)^2a ' 



(1) 



with the DGP model described by a = 1/2. More gen- 
erally, a is a freely-specifiable parameter, modulo two 
constraints: Firstly, in order for the modification to be 
relevant in the infrared, we must demand that a < 1. 
Secondly, it is straightforward to show that this propa- 
gator has a well-defined spectral representation, and is 
therefore ghost-free, only if a is positive definite [38]. 
Thus we have the allowed range: 



< a < 1. 



(2) 



The above parametrization is not only phenomeno- 
logically simple, but it also makes contacts with DGP 
(a = 1/2), as well as its higher-codimension generaliza- 
tions [31, 32, 33, 34]. Since the bulk has D>% space-time 
dimensions in the latter class of theories, the force law 
scales as l/r^~^ in the infrared, and the corresponding 
propagator tends to a constant {D > 6) or behaves as 
logfc (£> = 6) as fc ^ 0. Thus all higher-dimensional 
extensions of DGP correspond to a « theories [34]. 



A. Strong coupling of the longitudinal mode 

The propagator (1) describes a resonance, or a contin- 
uum of massive states, which is peaked at the tiny scale 
r~^. An immediate consequence is that the graviton now 
propagates 5 degrees of freedom: the 2 helicity-2 states 
of Einstein gravity, 2 helicity-1 states, and 1 helicity-0 or 
longitudinal mode, usually denoted by n. 

On top of the massive spin-2 representation, higher- 
codimension extensions of DGP also have D — 5 scalar 
fields that couple to T^. For simplicity, we will assume 
that these scalars have the same propagator as that of the 
massive graviton (1), although this is generally not the 
case in consistent realizations [33, 40]. In other words, for 
our purposes we will not distinguish between the helicity- 
mode and the extra scalars, and will denote them col- 
lectively as TT. 

These extra scalar degrees of freedom contribute to the 
trace part of the one-particle exchange amplitude [34]: 



1 



fc2 + ^-2(l-a)^2o 



rpfj.l/rp/ 



D-2 



-TT' 



(3) 



Because the resulting \/{D — 2) coefficient differs from 
the usual 1/2 of General Relativity (GR), this linearized 
amplitude is apparently inconsistent with solar system 
tests. And, remarkably, this departure from GR persists 
even in the massless limit, Tc ^ oo, a phenomenon fa- 
mously known as the vDVZ discontinuity [16]. 

As conjectured by Vainshtein in the context of mas- 
sive gravity [17], however, these theories can be phe- 
nomenologically viable because the linearized approxi- 
mation for TT breaks down in the vicinity of astrophysical 
sources [41, 42, 43, 44]. Below a macroscopic scale r*, 
non-linear interactions in tt become important and result 
in its decoupling. Power counting arguments [38] give 



l/(l-F4(l-a)) 



(4) 



where rsch is the Schwarzschild radius of the source. And 
since —> c» as rc —> oo, GR is recovered everywhere 
in this limit, as it should. Approximate explicit solu- 
tions [41, 45, 46] in DGP have confirmed the Vainshtein 
effect and the recovery of GR for r <i; r*. 

The leading correction to GR near a source can be 
inferred as follows. First of all, the solution for 7r(r) must 
be of order rgch/''* at r = r*, in order to match the linear 
solution. Furthermore, within r^, we expect the solution 
to be analytic in Tc These requirements together 

fix the parametric dependence of tt for r <C [38, 39] : 



^•Sch 



-H2(l-a) 



(5) 



In other words, the correction to the Newtonian potential 
is of order 



5^ 



r I r 



2(l-a) 



(6) 
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An important constraint on such modification comes 
from Lunar Laser Ranging experiments [39, 47]: 6"^/"^ < 
2.4 X 10"" [48]. Substituting rgch = 0.886 cm for the 
Earth and r = 3.84 x 10^° cm for the Earth-Moon dis- 
tance gives 



rc ^4-10" 



(7) 



This constraint on {a,rc) is shown in Fig. I. ExpHcitly, 
for standard DGP (a = 1/2) this gives > 120 Mpc. For 
a = 0, on the other hand, the bound is far weaker: rc ^ 
1 pc. Future lunar precession experiments will improve 
these bound by an order of magnitude [48]. Somewhat 
weaker constraints on DGP have also been obtained by 
studying the effect on planetary orbits [49], and it will 
be very interesting to extend these studies to general a 
theories. 



III. BACKGROUND COSMOLOGY 

The standard DGP model has a modified Friedmann 
equation of the form [50] 



SttG H 
3 rc 



(8) 



where we have assumed spatial flatness. Much atten- 
tion has been paid to the "plus" (or self-accelerating) 
branch of this equation, because of its asymptotically 
de Sitter solution even in the absence of vacuum en- 
ergy. However, various arguments have by now estab- 
lished that this branch suffers from ghost-like instabili- 
ties [43, 44, 51, 52, 53, 54]. We shall henceforth focus 
on the "minus" or healthy branch. Since the 1 /Hrc cor- 
rections slows down the expansion rate in this case, we 
must include a cosmological constant to account for cos- 
mic acceleration. 

A simple phenomenological extension of this equation 
suggests itself: 8nGp/3 - H"^-'' /rj. (The "plus" 

branch version of this equation was introduced in [55] 
to study generalized self-accelerated solutions.) We can 
fix 7 in terms of a as follows. First, from (8) we see 
that 7 = 1 for standard DGP, corresponding to a = 1/2. 
Furthermore, note from (1) that a = 1 can be absorbed 
into a redefinition of G\ if this is also the case for the 
Friedmann equation, then 7 = for a = 1. Assuming 
for simplicity a linear dependence of 7 on a, we conclude 
that 7 = 2(1 — a), or 



SttG , 



,2(l-a) 



(9) 



where A is a cosmological term. In the higher- 
dimensional framework, A should flnd a natural origin, 
for instance as a remnant vacuum energy component 
from the degravitation mechanism [39]. 

Remarkably, this parametrization implies that a = 
(corresponding to 2 or more extra dimensions) leads to 



identical expansion history as ACDM cosmology. Of 
course, higher-codimension DGP models will lead to 
small deviations from ACDM expansion history, as they 
do not exactly correspond to a = 0. But a = should 
offer a realistic approximation to their modified Fried- 
mann equation since we expect the correction to be a 
slowly- varying function of H/Vc [56] . 

For small values of a, we can quantify the expansion 
history in terms of the usual dark energy parameters 



W,'Wa 



PDE.eff 



87rG 



A- 



2(l-a) 



W 



-1 
-1 



1 dlnpDE,eff 
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din a 

{Hrc)Hl-n^) ' 



din a {Hrc)'^{l-n^) 



(10) 



Note that since a is positive, the modified gravity cor- 
rection leads to an effective dark energy component with 
w < —I. This is indeed the case in the DGP model [57]. 
While such phantom behavior may at first seem surpris- 
ing, it is in fact a natural consequence of the scalar-tensor 
nature of the effective theory in DGP. Indeed, at the 4D 
effective level the Friedmann equation is cast in Jordan 
frame, which can lead to an effective w < —1 [58, 59, 60]. 

The current (mainly) geometrical constraint w = 
—0.972 ± 0.060 (1(7 error), from a combination of 
WMAP5, baryonic acoustic oscillations, and super- 
novae [70], implies 1 + w > —0.067 at 95% confidence 
level if we impose the prior that w < —1, which leads to: 



> 3.2 [1 + ^("^ 



(11) 



In other words, unless a < 0.1, Vc should be larger than 
today's Hubble radius. 

Figure 1 compares the solar system constraint on r^, — 
see (7) — with the cosmological constraint, if we inter- 
pret (11) as a constraint on the distance to the last scat- 
tering surface (assuming a spatially flat universe). We 
see that for a < 0.5, the main constraint on the value of 
Vc comes from the background cosmology. However, as 
we noted above, in the a — > limit (corresponding to 2 
or more extra dimensions) , the background cosmology is 
indistinguishable from ACDM, and thus rc remains un- 
constrained. We will see in the next section that cosmo- 
logical perturbations will be sensitive to rc even in this 
limit. 



IV. COSMOLOGICAL PERTURBATIONS 

Since the study of cosmological perturbations in ex- 
plicit higher-codimension extensions of DGP has yet to 
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a 



zu = —2g/{l + g). A key difference is that these pa- 
rameters were all assumed to depend on a only, whereas 
our g(a, k) will also be a function of scale. 

For future reference, it is convenient to introduce the 
lensing potential 



Gravitational lensing observations, as well as the Inte- 
grated Sachs-Wolfe (ISW) effect, are both sensitive to 
Galaxy peculiar velocity measurements, on the 
other hand, are determined by the Newtonian potential 
\E' = (g — 1)(1>_. See [61] for a pedagogical description of 
how various observables relate to different metric combi- 
nations. 

In the next few subsections, we will specify g(a, k) for 
our modified gravity theories, based on known results 
in standard DGP. An important lesson from studies of 
DGP is that g{a, k) has qualitatively different behaviors 
on sub- or super-horizon scales. Hence we will explore 
these two regimes separately. 



FIG. 1: 95% lower limits on Tc, the scale of transition to higher 
dimensional gravity, from background cosmology (see (11)) 
and lunar laser ranging (see (7)), as a function of a. Standard 
DGP corresponds to a = 1/2, whereas its higher-dimensional 
extensions all have a w 0. 



be worked out, for our analysis we must resort to a phe- 

nomenological parametrization of these modifications. In 
this paper we take the parametrized post-Friedmann 
(PPF) approach introduced recently [37]. This formal- 
ism is reviewed in Appendix A. 

In Newtonian gauge, scalar metric perturbations are 
specified by the gauge-invariant potentials ^ and 



-(1 + 2*)dr + a\l + 2^)dx 



(12) 



In the absence of anisotropic stress, the momentum con- 
straint in GR fixes $ = However, this is no longer 
true in DGP because of the helicity-0 mode. At the 4D 
effective level, this can be understood tracing back to 
the scalar-tensor nature of the theory, as it is well-known 
that $ ^ in Jordan frame. This deviation from GR is 
parametrized in PPF through 



g{a,k) 



$ -I- * 



(13) 



This function can be specified independently from the 

expansion history. 

At this point we should mention closely related 
parametrizations of modified gravity that have appeared 
in the literature. Jain and Zhang [61] parametrized the 
deviation from $ = — \E' through 77 = 1 + g/{l — g). 
This convention was also adopted by Bertschinger and 
Zukin [62], except that the parameter was denoted by 
7. Daniel et al. [63], on the other hand, introduced 



Sub-horizon evolution 



Focusing on sub-Hubble scales and non-relativistic 
sources in standard DGP, Lue et al. [22] obtained 



5DGP|fe»aH 



3 [l + 2ifre(l+^)] 



(15) 



where ' = d/dlna. Note the plus sign within the brack- 
ets, which is appropriate for the normal branch [29]. 
Apart from this sign, this result was also obtained in 
cosmological perturbation theory using a quasi-static ap- 
proximation [24]. 

The dependence on Hrc in (15) can be attributed to 
the dynamics of tt. Indeed, for the universe as a whole, 
both r and rsch can be approximated by H~^. Setting 
a = 1/2 in this case, (6) gives the leading correction to 
standard cosmology, (5^'/^'|jjQp ^ l/Hvc, in agreement 
with (15) in the strong coupling regime Hrc ^ 1- Re- 
peating the argument for general a, (6) gives 



(5* 



(ifrc)2(i-«) 



(16) 



Meanwhile, in the weak coupling regime, from the ex- 
change amplitude (3) we can read off (5\E'/\E' = —(£) — 
4)/(_D — 2), consistent with the numerical coefficient 
of (15) for D = 5. Taking into account this numerical 
coefficient, this suggests generalizing (15) to 



D-A 
D-2 



2 {Hr, 



\2(l-a) 



(1 



H' 

3H ^ 



(17) 

A similar parametrization of g for general a theories was 
also studied in [64], albeit for self-accelerated cosmolo- 
gies. Our parametrization is also related to the slip pa- 
rameter zu introduced in [63]. The precise translation, as 
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mentioned earlier, is m = — 2g/(l + g). In the matter- 
dominated era, these authors assumed vu ~ ^/H^, which 
coincides with the scaUng of g in the case a = 0. 



B. Super-horizon evolution and decay of 
Newtonian potential 

In standard DGP, the evolution of super-horizon fluc- 
tuations was studied in [25] using a scahng ansatz for 
bulk perturbations. On the self-accelerating branch, the 
solution is well-fitted by [37] 



(18) 



The above parametric dependence and the form of our 
sub-Hubble parametrization (17) suggest the following 
ansatz for general a, 



9\k<^aH 



l + 2{Hr.f^'-^ (1+^) 



(19) 



The normalization constant A will be fixed shortly. 

Using a simple interpolation between the self- 
accelerated analogue of the short-wavelength fifoGP given 
in (15) and its long-wavelength counterpart (18), Wang 
et al. [30] recently showed that the self-accelerated DGP 
model is disfavored at the ~ 5a level compared to ACDM. 
Their analysis includes CMB temperature anisotropy, su- 
pernovae data, and constraints on the Hubble constant. 
Roughly 70% of the tension comes from the modification 
in the Friedmann equation, while the remaining 30% is 
due to the difi^erence in growth history, resulting in too 
large an ISW contribution to the low-^ multipoles of the 
CMB. 

Similarly, we find that an unacceptably large ISW 
component is a generic outcome for general a theories, 

unless Tc is much larger than the present Hubble ra- 
dius. But there is one important and very interesting 
exception: by choosing A Ki 1 and making a few mi- 
nor assumptions on the remaining PPF parameters to 
be discussed shortly, the Integrated Sachs- Wolfe effect 



resulting from the modified growth history can cancel 
the primordial Sachs- Wolfe contribution, resulting in a 
power deficit on large angular scales. This novel mech- 
anism can therefore naturally account for the observed 
low quadrupole as well as the lack of angular correlations 
above - 60° [11, 12, 13, 14, 15]. See Fig. 3. 

The actual value of A in explicit higher-codimension 
modifications of gravity remains at present unknown and, 
in any case, is likcJy to be model-dependent. Thus we 
view the agreement with large-angular CMB anomalies 
for ^ = 1 as a strong observational guide towards build- 
ing successful models of infrared modified gravity. The 
onus is then on model-builders to find explicit brane- 
world constructions where A = 1 is realized. For present 
purposes, we shall assume this can be achieved and ex- 
plore the consequences for other observables. Let us 
therefore take 



9\k<s:aH 



(20) 



Note that the transition g 1 has a natural physical 
interpretation: from the definition of g in (13), this cor- 
responds to a decay of the Newtonian potential \E' relative 
to the lensing potential 

At this point the astute reader might worry about the 
sign oi g\k^aH- On sub-Hubble scales, a positive g would 
be interpreted as tt mediating a repulsive force and there- 
fore representing a ghost. Fortunately, such reasoning 
does not apply on super-horizon scales, and it is straight- 
forward to devise toy scalar-tensor theories in 3+1 dimen- 
sions that lead to a positive g on super-Hubble scales. For 
instance, as shown in Appendix B, all Brans-Dicke the- 
ories satisfy this property during the matter-dominated 
era. 



C. Summary of fiducial model 

The full g{a, k) must interpolate between the sub- and 

super-Hubble behavior given in (17) and (20), respec- 
tively. For this purpose, we choose the interpolation form 



g{a,k) 



2 (iJrc) 



2(l-a) 



D 



-9 aH 



.b -\ 



(21) 



where a, b and Cg are constants. While our results are not 
overly sensitive to the precise values for these constant 
parameters, in practice we have found that smoother in- 
terpolations offer a better fit to the low multipoles of the 
CMB power spectrum. Hence, as fiducial values we take 



a = 6 = 0.5 ; Cg = 0.5 . (22) 

Although the interpolation is smooth, one can think of 

Cg as setting the scale below which g becomes negative. 
Note that for the above choice of Cg, this scale is compa- 
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rable to the horizon. Larger values of a, & and Cg diminish 
the effectiveness of the cancelation mechanism. 

Another key input parameter in PPF is cp, which per- 
tains to the evolution of the curvature perturbation C 
(see (A2) for its definition). Independent of the theory 
of gravity under consideration, energy-momentum con- 
servation alone imposes that C, is conserved in the infinite- 
wavelength limit [65]. Using this fact, Bertschinger [66] 
derived the following consistency relation 

C = const. = ^ [{g - l)<i>_ - 5'*- - (5 + I)*'-] 

+(.9 + !)$_. (23) 

In general, this condition must hold well outside the hori- 
zon. As we will see shortly, however, here we impose that 
Q is conserved on all scales in our fiducial modified gravity 
model. 

While C is actually conserved on all scales in ACDM 
cosmology, this is not necessarily the case for modified 
gravity theories. In the PPF approach, one specifies a 
parameter cr which, as described in Appendix A, de- 
termines an effective horizon cr/aH above which ^ is 
conserved. Our mechanism for canceling the large-scale 
CMB power requires that conservation of C persist on a 
sufficiently wide range of sub-Hubble modes, correspond- 
ing to small values of cr- In practice, we have found that 
Cr ^ 10~^ is desirable. For simplicity, however, as 
fiducial value we shall set 

cr = , (24) 

thereby enforcing C conservation on all scales. The evo- 
lution of super-horizon modes is therefore uniquely de- 
termined through (23) by specifying g and the expan- 
sion history. Note that the assumption of constant ( was 
also made by Bertschinger and Zukin [62] in what they 
dubbed scale-independent modified gravity models. (A 
key difference, however, is that their g was also indepen- 
dent of scale.) In Appendix A, we shall describe the PPF 
formalism in more generality, allowing for non-zero cr- 

In light of the relatively large number of input func- 
tions and parameters inherent to the PPF framework, 
let us synthesize the essential qualitative and quantita- 
tive features of our fiducial parametrization: 

• The difference in metric potentials g{a, k) is cho- 
sen to interpolate between the sub-horizon behav- 
ior, dictated by well-understood features of mas- 
sive/resonance gravity outlined in Sec. II A, and 
the more model-dependent super-horizon behavior. 
For the latter, the large-scale CMB anisotropy con- 
strains us to choose A = 1, so that g\k^aH makes 
a transition from at early times to 1 at late times 
— see Fig. 2. Physically, this corresponds to a 
decay of the Newtonian potential ^' relative to <&. 
Smoother interpolations are preferable. 

• The curvature perturbation must be nearly con- 
served, which requires small values for cr- For sim- 
plicity, we set Cr = in our fiducial model. 
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FIG. 2: Dependence of g{k,a) (see (21)) in our model on 
k/aH in the matter-dominated era for = 0.3, 1,3, and 10 
Hubble radius. The number of space-time dimensions, D, is 
assumed to be 6 (i.e. 2 extra dimensions). General Relativity 
is recovered when rc ^-ff~^ as, g—>0. 

The remaining parameters to vary are of course a, D 
and re- Since we are interested in higher dimensional 
generalizations of DGP, we shall fix 

a = , (25) 

as discussed below (2). Since the expansion history re- 
duces to ACDM in this case, we are free to consider val- 
ues of Vc that are relatively small compared to the Hubble 
radius. We will find that our modified gravity theories 
can explain the anomalies listed in the Introduction and 
match other observations with Vc ^ 300 — 600 Mpc. Fi- 
nally, the total number of space-time dimensions D sets 
the strength of the extra scalar force at short distances, 
as seen in (17), with larger D corresponding to stronger 
gravity. To explain the CBI excess and the anomalous 
Lyman-a power, it will suffice to consider two extra di- 
mensions, or 

D^6. (26) 
This completes the description of our fiducial model. 

V. CMB TEMPERATURE ANISOTROPY AND 
LARGE-SCALE ANOMALIES 

Let us step back and explain why the choice for g\k<^aH 
in (20) can explain the low quadrupole and the lack of 
correlation on > 60° scales. On large angles, the CMB 
anisotropy is well approximated by the sum of the Sachs- 
Wolfe and the Integrated Sachs- Wolfe (ISW) effects [67]: 

= V+2/dt^^ V + 2AcI>_. (27) 
TcMB 3 J dt 3 ' 

The integral in the ISW effect is taken along the light 
cone, while the partial derivative is at fixed comoving spa- 
tial position. Therefore, the last approximation in (27) 
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is only valid when the spatial variations of are much 
smaZZer than its time variation: |V<1>_| <C It is then 

interesting to note that on these scales, A<I>_ can be ob- 
tained analytically as a function of g(a, k), from the fact 
that the curvature perturbation ^ remains constant on 
super-horizon scales — see (23). In particular, if g goes 
from to a finite (constant) value during the matter 
era, from (23) we see that: 



9i 



(28) 



For example, for the transition from to = 1 we have: 



STc 



MB 



1 



Tc 



2A$_ = 



MB 



5(1 - 9i) 



3(5 + <?f) 



$_=0! (29) 



In other words, if (/ — > 1 at late times in the modified 
gravity theory, the ISW and Sachs- Wolfe effects exactly 
cancel on large angles — corresponding to super-horizon 
scales at the time of the transition — , which could po- 
tentially explain the vanishing of CMB power on scales 
> 60°. Note that g ^ 1 is equivalent to a vanish- 
ing Newtonian potential relative to the lensing potential, 
'i' <C at late times. A partial cancelation between 
the primordial Sachs- Wolfe and a late-time ISW contri- 
bution arising from g was also noticed in [62, 63]. 

This phenomenological observation, in addition to the 
fact that Brans-Dicke theories have g > on super- 
horizon scales — see Appendix B — , justifies the infrared 
limit of our ansatz for g(a, k) given in (21). 

Figure 3 demonstrates how this cancelation can work 
at low Cs for different values of r^. The CMB power 
spectra were generated using a version of CAMB [68] 
modified by W. Fang to include PPF parameters [69]. 

The apparent lack of power on large angles was first 
quantified by Spergel et al. 2003 in the first data release 
of WMAP [71], through the quantity: 



r.l/2 

/ [C(0)]2dcos0, 



(30) 



where C{9) is the angular correlation function of the CMB 
temperature anisotropies shown in Fig. 4. Most recently, 
Copi et al. [13, 15] argued that the probability of S1/2 to 
be lower than it is observed outside the Galactic plane is 
only 0.025% in the concordance ACDM model. However, 
different estimators for the correlation function can lead 
to larger values of 6*1/2. Figure 5 compares this proba- 
bility for our modified gravity models with ACDM, and 
shows that it is increased by about a factor of 2-3 in 
the relevant range. While this is not enough to explain 
the low values of 5*1/2 in direct measurements of corre- 
lation function [e.g., the two smaller values in Fig. 5], 
it can bring the prediction within 90% probability range 
for maximum likelihood estimates of the WMAP5 power 
spectrum [the larger value in Fig. 5]. 

Despite the improvement of the consistency be- 
tween models and data at large angles for temperature 
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FIG. 3: CMB angular power spectra for best-fit ACDM 
(blue, solid curve), — 600 Mpc (green, dashed curve) and 
rc = 300 Mpc (red, short-dashed curve). The vertical bars 
show the cosmic variance spread centered on the respective 
theoretical predictions. In the middle and bottom panels, the 
previous curves are shown for comparison. The data points 
are from the WMAP-5 data release [70] . 



anisotropy, the marginally increases by 1.21 (0.32) for 
Tc = 300 (600) Mpc, with other cosmological parameters 
fixed. The origin of this increase can be seen in Fig. 6, 
where the suppression of power at small ^'s decreases 
the TE temperature-polarization cross-power spectrum 
below the observed value. However, given that the low- 
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FIG. 4: Square of the correlation function of CMB tempera- 
ture anisotropy, as a function of the cosine of the separation 
angle in the sky. The curves show ACDM (solid), = 300 
Mpc (dotted) and = 600 Mpc (short-dashed). The long- 
dashed curve is the Legendre transform of the WMAP5 max- 
imum likelihood power spectrum [72]. The observed correla- 
tion is systematically below the ACDM prediction for > 60° , 
or 0036* < 0.5. 
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FIG. 6: Temperature-Polarization (TE) and polarization 
(EE) power spectra for ACDM (blue, solid curve), = 300 
Mpc (red, short-dashed curve) and — 600 Mpc (green, 
dashed curve). The points are from the WMAP5 data re- 
lease [72]. 



est ^-bin in the TE spectrum decreased by a factor of 
3 from the 1st year to the 5-year WMAP data release 
due to an improved foreground cleaning, it is still con- 
ceivable that the systematics in this bin are underesti- 
mated. Therefore, we predict a significantly lower TE 
cross-power spectrum at £ < 10, which should be clearly 
distinguished from ACDM by the Planck satellite, due to 
its better polarization sensitivity and foreground cleaning 
capabilities [73]. 

Meanwhile, as shown in the bottom panel of Fig. 6, the 
EE spectrum is not appreciably affected by our modifi- 
cation of gravity. 



FIG. 5: The cumulative probability for 5*1/2, defined in (30), 
measuring the CMB correlation on scales > 60°. The curves 
are for ACDM (solid), = 300 Mpc (dotted) and = 600 
Mpc (dashed). The vertical lines show the observed values of 
Si/2 from different estimators (see [15] for details). 



We should also point out that it is possible to improve 
the full (relative to ACDM) without extra marginal- 
ization, simply by decreasing b in (21). However, this 
will spoil our prediction of extra growth for small-scale 
structures, which we shall discuss next. 
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FIG. 7: DifTerent growth regimes as function of scale factor 
and comoving wavenumber k. Colored regions correspond 
to modes inside the horizon today. The solid line labeled 
by aH denotes the Hubble horizon; the dotted line traces 
the characteristic scale set by Cg. At early times, H > r~^, 
growth proceeds as in GR (Green region); for H < r^^ and 
on scales much smaller than the Cg-scale, growth is enhanced 
thanks to the helicity-0 mode (Blue region); for H < 
and on large scales, the decay of the Newtonian potential ^ 
suppresses growth (Red region). 



VI. STRUCTURE FORMATION WITH 
MASSIVE GRAVITY 

The growth of perturbations in our model can be un- 
derstood analytically by making a few simplifying ap- 
proximations. The essential physics is illustrated in 
Fig. 7. The solid line labeled by aH denotes the Hubble 
horizon. Meanwhile, the dotted line denotes the scale set 
by Cg where g makes a transition from positive values on 
larger scales to negative values on smaller scales. 

At early times, H ^ , we have a GR regime on 
all scales, as « 0. For modes much smaller than the 
Cg-scale, once H drops below the scalar force medi- 
ated by the longitudinal mode kicks in, enhancing the 
growth of perturbations. This is the blue-shaded region 
in Fig. (7). Thus there is excess of power on small scales, 
compared to ACDM. 

For modes with wavelength above the Cg-scale, the 
function g makes a transition from to 1 at _ff '--^ r~^, 
corresponding to a decay of the Newtonian potential 
This is the red-shaded region in Fig. 7. This suppression 
of growth results in a power deficit in the spectrum of 
density perturbations compared to ACDM. 

Before numerically solving for the evolution of pertur- 
bations, we can gain analytic insights by studying dif- 
ferent regimes. Although our perturbed metric is cast 
in Newtonian gauge, we find it more convenient to track 
comoving gauge density perturbation, A^. Its evolution 
equation follows as usual by combining energy and mo- 
mentum conservation [see (Al)] 

A;;+(2+f )A;, = (l-5)^<i>», (31) 

where we have used the fact that the curvature perturbar- 
ion is assumed to be conserved on all scales, C,' = 0. As it 



should, (31) reduces to the standard evolution equation 
in the GR limit Tc — s- oo. 

At late times {H < ?'^^), when the transition is com- 
plete, g{k^ a) only varies slowly as the mode wavelength 
relative to the horizon shrinks with time — see Fig. 2. 
In this regime we can ignore <!>'_ and g' in the C, equa- 
tion (23), which leads to (28) in the matter era. With 
this in mind, we can view (31) as a forced damped har- 
monic oscillator, where the force is enhanced by a factor 
of 5(1 — g)/{b + g) compared to ACDM. Therefore, prior 
to the domination of A the inhomogeneous solution be- 
comes 

^_^^^CUM. A^^^il^A^DM. (32) 

b+g b+g 

On very small scales, where g ^ —{D~A)/{D~2), this 
gives 

j>_ ^ 5p-2) j>CDM 5(i^-3) ^cDM^ 

2(215-3) " 2I?-3 " 

(33) 

In particular, for our fiducial case of D = 6, 

^ ^ci,CDM . A„,^^ACDM^ (34) 

Meanwhile, on large scales we have g — s- 1, and thus 
Am ^ while <i>_ decreases by a factor of 5/6, as we 
discussed in the previous section. 

Figure 8 demonstrates this behavior for and Am 
with D — 6 over the range k = lO"'* — 10 Mpc^^ cover- 
ing both super-horizon and sub-horizon modes. (What 
is plotted in the right panel can be thought of as the 
transfer function for an effective potential related to Am 
through the standard Poisson equation.) As expected, 
prior to the domination of A, large scale (super-horizon) 
modes decay, while small-scale sub-horizon modes grow 
with time. 

Once A dominates, however, all modes start to decay. 
Indeed, since ff' — > in this regime, from (23) we find: 

-I- (^-[^^ = =^ cx a^Tff . (35) 

Thus, while <i>_ decays as a~^ at late times in ACDM 
cosmology, here it decays even faster, as a~^^^, since 
g — —{D — A) / {D — 2) on small scales in our model. In 
particular, for D = 6, we have <!>_ ^ a~^ at late times. 
This behavior can be seen in Fig. 8, where initially 
grows by up to ^ 10% as the graviton Compton wave- 
length enters the horizon, in agreement with (34), but 
then decays faster than in the ACDM model after the 
onset of cosmic acceleration. The competition between 
these two effects is why, as we will see below, gravita- 
tional weak lensing is not substantially affected by the 
graviton Compton wavelength in our model at low red- 
shifts. 

In the following sections, we will discuss the predictions 
of our model for different measures of structure formation 
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FIG. 8: Transfer functions for the lensing potential <!>_ and the comoving density perturbations Am, as a function of redshift 
for ACDM (sohd), n = 300 Mpc (dotted) and rc = 600 Mpc (dashed), all assuming D = 6. The curves show k = 10""*, 3.6 x 
10"*, 1.3 X 10"^, 4.6 X 10"^, 1.7 X 10"^, 6 x 10"^, 0.2, 0.8, 3, and 10 Mpc"^ from bottom to top, and are displaced for clarity. 
Prior to the domination of A, large scale (super-horizon) modes decay, while small-scale (sub-horizon) modes grow with time. 
Once A dominates, all modes start to decay with time. 
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FIG. 9: The power spectrum of the linear comoving density 
Am at 2: = 0.1 (solid), z = 1 (dotted), and z = 5 (dashed), for 
Tc = 300 Mpc. For clarity, the spectra are normalized to unity 
at fc = 10 Mpc~^. We see that, unlike the ACDM model, the 
shape of the power spectrum evolves on large scales, which 
will cause the galaxy bias to be generically scale-dependent. 



at low redshifts. One notable omission is the clustering 

of galaxy redshift surveys, such as SDSS and 2dF, which 
have been instrumental in measuring the shape of the 
matter power spectrum at low redshifts. These studies 
T(Ay on a key assumption: that galaxy distribution follows 
matter distribution on large scales, up to an arbitrary 
but constant bias fac;tor. Unfortunately the assumption 
of constant bias only holds if the shape of the linear mat- 
ter power spectrum does not change with time, as is the 
case in the ACDM model. Here, on the other hand, the 
shape of the matter power spectrum does evolve in time, 
as shown in Fig. 9. As pointed out in [74], we there- 
fore expect a scale-dependent bias, since the shape of the 
power spectrum of galaxies may depend on that of mat- 
ter during the entire history of galaxy formation. We 
shall postpone a full study of this effect, and only warn 
that even the large scale bias in modified gravity models 
cannot be captured by a single number. 



A. Weak lensing power spectrum 

Weak lensing is one of the few direct probes of the 
large-scale gravitational potential with minimal confu- 
sion from complicated astrophysical processes. While 
early weak lensing measurements were plagued by ob- 
servational systematics, there have been much recent im- 
provement in controlling these errors. Here we will fo- 
cus on the recent measurement of large-angle weak lens- 
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FIG. 10: The ratio of lensing correlation measured within 100' 
radius for Vc = 300 Mpc (dotted) and Vc = 600 Mpc (dashed) 
with £) = 6, to the ACDM prediction. The data point is 
from CFHTLS Wide weak lensing measurements [75], with 
the vertical error bar denoting the Icr range. The horizontal 
error bar spans the redshift range for the relevant lensing 
kernel. 



ing correlations with Canadian Prance Hawaii Telescope 
(CFHT) observations [75], which finds: 

asf^j = 0.837 ±0.084 (36) 

on angular scales 85' < 9 < 230'. Comparing this 
with the WMAP5 combined constraints with baryonic 
acoustic oscillations (BAO) in galaxy surveys and type 
la supernovae [70], we see that the lensing amplitude is 
0.98 ± 0.11 times the CMB-preferred value, i.e. they 
are in excellent agreement. Figure 10 compares this re- 
sult with predictions of our massive gravity model. To 
model the CFHT weak lensing measurements, we plug 
in our model into the standard Poisson equation, 
and compute the resulting density fluctuations within a 
sphere of angular radius 100', roughly corresponding to 
the scale of weak-lensing measurements. As the CFHT 
lensing sources are broadly distributed in the redshift 
range 0.5 — 2, we assume the range z = 0.25 — 1 for 
the lensing kernel. 

Therefore, due to the competition between massive 
gravity effects and cosmic acceleration discussed above, 
we see that current weak lensing measurements can- 
not distinguish our predictions from that of the ACDM 
model. However, future percent-level tomographic mea- 
surements of weak lensing power from the Large Synop- 
tic Survey Telescope [76], Pan-STARRS [77], or NASA's 



FIG. 11: Cross-power spectrum of a galaxy survey at z ~ 0.5 
with the CMB, for ACDM (soUd), Vc = 300 Mpc (dotted) and 
Tc = 600 Mpc (dashed), again for D = 6. The data point is a 
rough measure of the current observations [2] (see the text). 
The horizontal error bar spans the relevant angular range. 

Joint Dark Energy Mission [78] experiment will be able 
to test our predictions for the lensing potential. 



B. Integrated Sachs-Wolfe Cross Correlation 

The cross-correlation between the ISW effect in the 

CMB and the large scale structure at low redshifts probed 
by galaxy surveys, was first proposed in [79] and later de- 
tected by several groups in cross-correlation with WMAP 
CMB maps. Most recently, [2] and [3] combined different 
galaxy surveys to obtain an overall detection of the ISW 
effect in cross-correlation at the 4 — 4.5f7 level, within 
the redshift range z = 0.1 — 1.5. However, it was also 
pointed out in [2] that the amplitude of this detection 
is 2.23 ± 0.60 larger than what is predicted by ACDM 
model (see also [80] for a similar result, but at lower sig- 
nificance). Going back to Fig. 8, we note that a faster 
decay of $_ is indeed expec;ted in our mock^l at low red- 
shifts. In Fig. 11, we compare our prediction for the ISW 
cross-power spectrum for a galaxy survey at ^; ~ 0.5, with 
that of ACDM. The data point, associated with the ex- 
cess observed amplitude, is centered at £ ~ 20, as most 
of the S/N in ISW detection in cross-correlation comes 
for z - 0.5 and £ - 10 - 40 [81]. Since the ISW effect is 
proportional to as seen in (27), we have scaled the 
ACDM prediction for cross-power by the ratio of for 
our model to the standard one, in order to obtain the 
dotted and dashed curves in Fig. 11. Additional correc- 



12 



\ 2 




WMAP5+ 

0.8 



1.0 



1.2 



1.4 



1.6 



FIG. 12: The relative vs cts in ACDM and our modi- 
fied gravity models for producing the observed bulk flows 
on ~ 100 Mpc scales [4]. The vertical line shows the 
WMAP5+BA0+SN best fit value. 



tion due to the changes in the galaxy overdensity (which 
is correlated with the ISW) is possible, but since the bias 
is obtained by fitting the observed galaxy clustering in 
the present model comparisons [2, 3], these corrections 
are expected to be small. 

Interestingly, as can be seen from Fig. 11, the signal 
is fairly insensitive to the value of the graviton Compton 
wavelength over the range of interest. This is because, for 
these values of Tc, the transition to the modified gravity 



regime when H 



is already completed by z = 0.5. 



Therefore, we see that massive/extra-dimensional gravity 
can explain the observed excess in the ISW effect detec- 
tion in cross-correlation. 



C. Bulk Flows 

Another measure of structure growth on large scales 
are peculiar velocity measurements, which, through the 
continuity equation, probe in the linear regime. As 
discussed in the Introduction, recent measurements of 
coherent peculiar velocity (or bulk flows) on scales of 
'--^ 100 Mpc, from very different methods [4, 5], have led 
to consistent magnitudes and directions which are nev- 
ertheless significantly larger than the ACDM prediction. 
In particular, [4] find a mean velocity of 407 ± 81 km/s 
within a gaussian filter of radius 50 h~^Mpc, centered on 
the Milky Way, while the ACDM WMAP5 best-fit model 
predicts an r.m.s. value of 190 km/s. Figure 12 shows 
how the {— ^2 log of likelihood) for this measurement 
depends on the erg value of the ACDM model, keeping all 
other parameters fixed to their WMAP5 best-fit values. 
We see that the standard value — 0.812 ± 0.026 is 
disfavored at about 2(t. However, both of our fiducial 
models with Tc = 300 and 600 Mpc increase the r.m.s. 
velocity on 100 Mpc scale by 21%, which brings down 
the level of disagreement to 1.4a, as shown in Fig. 12. 
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FIG. 13: The transfer function of the density potential (re- 
lated to Am through the standard Poisson equation) at z = 3 
and k = 1.4 Mpc"^ for ACDM (sofid), = 300 Mpc (dotted) 
and Vc = 600 Mpc (dashed). The data point (with la error 
bar) characterizes the excess power observed in Lyman-a for- 
est observations [10]. The horizontal error bar approximates 
the range of scales probed by these observations. 



D. Lyman-a forest 

The absorption spectrum of the Lyman-a transition in 
the intergalactic neutral hydrogen on the line of sight to 
high redshift quasars, commonly known as the Lyman- 
a forest, has so far provided us with our most stringent 
constraints on the small-scale power spectrum of mat- 
ter perturbations at z — 2 — 3. The interpretation of 
these observations is complicated by mild non-linear evo- 
lution, as well as baryonic physics and other features in 
the quasar spectra, and thus large suits of numerical sim- 
ulations have been used to calibrate the relation of the 
observed flux spectrum to the underlying linear matter 
perturbations [8, 82]. In this section, we compare the 
predictions of our model to the constraints obtained on 
the linear theory power spectrum at z = 3 and fc ~ 1.4 
Mpc^i in [9]. Seljak, Slosar & McDonald [10] used these 
constraints to conclude that the predicted linear power 
spectrum from the WMAP3 best-flt model was too low to 
explain the observed spectrum of Lyman-a forest at the 
~ 2.4(7 level. They attributed this excess to the presence 
of extra neutrinos (or other relativistic components) dur- 
ing radiation-matter equality. Their best fit model with 
extra neutrinos has roughly 35 ± 14% more power than 
the standard ACDM WMAP3 best fit model with 3 neu- 
trinos. While a consistent interpretation of this discrep- 
ancy in the context of our model requires a similar suit of 
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FIG. 14: The ratio of as in our modified gravity models to 
that of ACDM, as a function of redshift. The dotted and 
dashed curves correspond to rc = 300 and 600 Mpc respec- 
tively. The data point shows the observed excess associated 
with the CBI excess [7]. 



simulations with a full treatment of non-linear modified 

gravity, as a crude first step Fig. 13 compares this exc;ess 
with our model prediction for the small-scale growth in 
linear density Am. Interestingly, either fiducial values of 
the graviton Compton wavelength in our model (rc = 300 
or 600 Mpc) can explain the observed excess. 



E. Galaxy clusters and CBI excess 

Finally, we will address another long-standing puzzle 
in CMB observations, known as the CBI excess. The 
high resolution observations of the Cosmic Background 
Imager (CBI) have long indicated an excess in the power 
spectrum of CMB anisotropies at ^ ~ 3000, far above 
the expected level of primary anisotropies [6]. Con- 
sistent results have been obtained by BIMA [83] and 
ACBAR [7] experiments, albeit at lower significance lev- 
els. While this excess has a consistent spectrum with that 
of the thermal SZ effect — the dominant source of CMB 
anisotropy at ^ ~ 3000 — , the observed level is signifi- 
cantly higher than what is expected in the WMAP5 best- 
fit ACDM cosmology. The amplitude of the SZ power 
spectrum is a direct measure of the abundance of mas- 
sive galaxy clusters at z ^ 0.5 — 2, and scales as erg [84]. 
The current excess in the SZ power can be interpreted as 
a 17% larger value of o-g at 2 — 3a level [7]. As ag is the 
amplitude of the linear matter density fluctuations close 
to the scale of cluster formation, we expect to capture 



most of the physics relevant for cluster abundances by 
tracking the evolution of (7g,(z) in our model. Figure 14 
shows our prediction for the excess of (Tg,(z) compared to 
the ACDM model, which could easily explain the CBI 
excess. 

We should note that the observed deficit of the SZ 
decrement in WMAP clusters [85] — indicating a smaller 
SZ fiux per cluster — may point to an even larger excess 
required to explain the CBI excess. Future observations 
of SZ power spectrum as well as SZ clusters by the At- 
acama Cosmology Telescope (ACT) and the South Pole 
Telescope (SPT) should be able to decisively confirm or 
rule out this excess. 

A seemingly contradictory result comes from the mass 
function of galaxy clusters at lower redshifts, inferred 
from X-ray observations. Most recently, Vikhlinin et al. 
[86] showed that Chandra observations of a sample of 
galaxy clusters detected in X-rays by ROSAT are remark- 
ably consistent with the predictions for mass function 
from ACDM, and in particular erg = 0.77 ± 0.02 (includ- 
ing systematic errors, but using geometrical constraints 
on flrn)- If correct, this would be inconsistent with the 
common interpretation of the CBI excess as a large erg 
[84], and may appear to rule out our preferred range of 
rc = 300 - 600 Mpc (Fig. 14). However, the relationship 
between the ROSAT X-ray luminosities and the cluster 
masses has a large scatter, which is not constrained very 
well empirically (see Fig. 13 in [87]). In particular, an 
extended low-luminosity tail in the distribution of lumi- 
nosities for a given mass would mean that a significant 
number of massive clusters may go undetected, leading 
to a possible underestimate of ag. 

We also should point out that our analysis so far is 
only strictly valid in the linear regime and the formation 
of non-linear collapsed objects (such as galaxy clusters) 
is not modeled in our framework. In other words, the 
relationship between the number of clusters and the lin- 
ear power spectrum, used in [86] , is only valid for ACDM 
cosmology, and could be significantly different in modi- 
fied gravity models. Indeed, as we discussed in II A, this 
is expected due to the Vainshtein effect at large densities, 
and is required in order to satisfy the precision tests of 
gravity at solar system and galactic scales. Therefore, 
a good understanding of non-linear collapse process and 
Vainshetin effect in our model is necessary to make reli- 
able predictions for the cluster mass function (and CBI 
excess) . 



VII. CONCLUSIONS 

Arguably the most pressing question in cosmology is 
whether the observed cosmic acceleration is due to a 

breakdown of Einstein gravity or to conventional vac- 
uum energy. Fortunately, the ACDM model is highly 
predictive — a single dark energy parameter A suffices to 
predict the entire expansion and growth histories. Near 
future probes of the large scale structure will subject this 
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standard model to increasingly stringent tests of its pre- 
dictions. 

In this paper we have presented a case that current 
observations may already be pointing towards new grav- 
itational physics on cosmological scales. The excess of 
power seen in the Lyman-a forest and small-scale CMB 
experiments, the larger bulk flows seen both in pecu- 
liar velocity surveys c md in kinetic SZ, and the higher 
ISW cross-correlation all present potential challenges for 
ACDM. And although their individual statistical signif- 
icance is only at the 2a level, it is certainly intrigu- 
ing that these observations suggest a common explana- 
tion: that structure is more evolved than expected from 
ACDM. On the other hand, weak Icnsing measurements 
give a value for ag consistent with the CMB. 

We have argued that these observations find a natural 
explanation in models where the graviton is a resonance 
with a tiny width Tc ~ 300 — 600 Mpc. The longitudinal 
mode TT of the graviton can be thought of as mediating an 
extra scalar force which, thanks to the Vainshtcin effect, 
only kicks in at late times and on sufiiciently large scales. 
This speeds up structure formation, thereby accounting 
for the above anomalies. Meanwhile, since photons do 
not couple directly to tt, weak lensing is much less af- 
fected, which explains its consistency with the CMB. 

Our approach has been motivated by higher- 
dimensional generalizations of the DGP scenario, such 
as Cascading Gravity models. A crucial advantage of 
considering more extra dimensions is that the modifica- 
tions to the Priedmann equation become indistinguish- 
able from a cosmological constant. This makes it possible 
to consider much smaller values of the graviton Compton 
wavelength than usually allowed by constraints on the ex- 
pansion history. And this in turn amplifies the impact of 
the longitudinal mode on structure formation. 

The obvious danger with enhancing structure growth 
is that one risks boosting the ISW contribution to the 
CMB low-£ power. In fact, observations show the oppo- 
site trend, starting with the well-known low quacirupolc. 
A more statistically significant anomaly is the peculiar 
fact that the angular correlation function vanishes on 
> 60° scales, as seen both in COBE and WMAP. 

We have shown that our modified gravity models can 
account for these large scale CMB anomalies. With cer- 
tain assumptions about the super-horizon behavior of tt, 
we have presented a novel mechanism in which a late-time 
ISW contribution cancels against the primordial Sachs- 
Wolfe amplitude. This ISW contribution is different from 
the usual one due to A, but arises because the dynamics 
of TT undergo a transition when H ^ r^^. Our mechanism 
relies on conservation of ( and is analogous to the jump 
in the Newtonian potential incurred during the radiation- 
to matter-dominated era in standard cosmology. 

Many issues remain to be addressed: 

• To compare our predictions with galaxy surveys, 
such as SDSS, we need a handle on the halo bias in 
our model. As argued in [74], the bias is generically 
expected to be scale dependent in modified gravity 



theories. This issue is currently under study using 
N-body simulations [88]. 

• A more thorough comparison with data requires a 
full likelihood analysis, including the PPF param- 
eters used here. 

• The effects of modified gravity is suppressed by the 
Vainshtcin effect at large densities. How does this 
affect the non-linear collapse process, and the halo 
mass function in our model? 

• While the phenomenological approach adopted 
here constitutes an important and instructive first 

step, it is imperative to justify the choice of PPF in- 
put functions through full-fledged brane-world cal- 
culations, for instance in Cascading Gravity mod- 
els. In particular, the suppression of the large-angle 
CMB power relies on certain assumptions about the 
behavior of perturbations on super-Hubble scales, 
as mentioned above, and it will be essential to check 
whether this behavior can be explicitly realized in 
higher-dimensional constructions. If not, are there 
other consistent theories of modified gravity that 
can reproduce the phenomenology presented here? 
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APPENDIX A: REVIEW OF PPF FORMALISM 

In this Appendix we review the essential ingredients 
of the PPF formalism, considering the more general 
case where ( is not necessarily conserved on sub-Hubble 
scales. Following [37] we study the matter perturba- 
tions in comoving gauge. In this gauge, the perturbed 
stress tensor for dark matter has non-zero components 
Tq = —p{l + Am) and = —pdiVm/ka, where Am and 
Win are by definition the energy density and velocity per- 
turbations, respectively. 

The matter stress-energy tensor is assumed to be co- 
variantly conserved, as is the case in the DGP model 
and its higher-codimension extensions. Conservation of 
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energy and momentum then implies 

k 



A' 



-3C'; 



(Al) 



with ' = d/dlna. Here ( denotes the curvature pertur- 
bation, related to the other variables by 



(A2) 



This can then be combined with (A4) to derive a second- 
order equation for 
3. Interpolating between sub- and super-horizon 
regimes 

In general, of course, the resulting sub-horizon evolu- 
tion equation for <i>_ will be inconsistent with conserva- 
tion of when applied in the long- wavelength limit. One 
must therefore interpolate between the sub- and super- 
horizon behaviors derived above. To do so, following [37] 
we introduce a new function r(a, k) through a modified 
Poisson equation, valid on all scales: 



1. Super-horizon evolution 



+r) =47rGaVAn 



(A7) 



Assuming adiabatic initial conditions, the curvature 
perturbation, must be conserved for long wavelengths. 

As shown by Bertschinger [66], this yields a consistency 
relation valid on super-Hubble scales: 

C = const. = §;[i9- 1)^- - q'^- -{9 + 1)^-] 



+(5 + 1)^- 



(A3) 



It follows that the evolution of super-horizon modes is 
uniquely determined by specifying g and the expansion 
history. 



At long wavelengths, k <C aH, we require that F 
enforces the conservation of By taking the time- 
derivative of (A7) and combining with the conservation 
equations (Al), we obtain 



r' + r = - 



Ml + /c) 



aH 



where specifies the leading behavior of (' as k ^ 0: 



1 



(A8) 



(A9) 



Sub-horizon evolution 



Then, taking the time-derivative of (A2), we can substi- 
tute for to obtain an equation valid for k <^ aH: 



On sub-Hubble scales, k aH. we enter the Newto- 
nian regime, in which metric and matter perturbations 
are related by a modified Poisson equation 



fc2 [1 + f{a, k)] = 47rGo^A„ 



(A4) 



Matching the above form factor with our modified prop- 
agator (1) fixes f{a,k): 



f{a,k) = 



kVf 



2(l-a) 



(A5) 



Note that / has been neglected in most earlier studies 
of DGP cosmological perturbations [22, 24, 29], except 
in [28] , presumably because Vc was assumed to be of or- 
der Hq^, and therefore all modes of interest satisfied 
k ^ In this work, however, we have considered 

much smaller values for rc, in which case /(a, k) must be 
carefully taken into account. 

On the matter side, the energy-momentum conserva- 
tion equations (Al) together yield an evolution equation 
for Am, valid in the Newtonian regime: 



H' 



a'^H'^ 



(l-ff)$_ 



(A6) 



V + T 



25 



1 
aH 



H' 



Meanwhile, it is clear from (A4) that F —^ f{a,k)^- 
at short wavelengths (fc ^ aH). To interpolate be- 
tween this condition and the long-wavelength behavior 
described by (AlO), we can assume that F satisfies 



1 + c; 



_aH 



(F-/(a,fc)$_) 



H' 



{l + g)H 2 



f^m(l + /c) • (All) 



Note that this is very similar to the corresponding equa- 
tion in [37], with the only difference that we are using 
{crk/aHY as opposed to (crk/aH)"^. A higher power of 
k is essential here to win over the potential l/k'^ depen- 
dence of /(a, k), achieved for a = 0, and thereby ensuring 
the correct infrared behavior (fc — > 0). 

Moreover, we can combine the conservation equa- 
tions (Al) and the Poisson equation (A7) to obtain 
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^(1 + 5)1 = ^(25-5')^- 



(r' + r) + ^ 



H' 



H k 



aH 



r 



(A12) 



Equations (All) and (A12), together with momentum 
conservation (Al), form a coupled system of first-order 
equations for F and v^. 

To summarize, the PPF formahsm requires specifying 
an expansion history H{a), a difference in metric poten- 
tials g{a,k), and a short- wavelength graviton form fac- 
tor f(a,k). Moreover, the transition between sub- and 
super-Hubble regimes is specified by cr, and the lead- 
ing departure from conservation of C is parametrized by 
In the particular case cp = = 0, which was con- 
sidered in the main text, is exactly conserved on all 
scales. Then, the evolution of is uniquely determined 
by (A3). Meanwhile, (A6), which holds exactly in that 
case, describes the evolution for A™. 



APPENDIX B: BRANS-DICKE 

To reassure ourselves that it is perfectly sensible to 
choose 5 > on super-horizon scales, as assumed in (20), 
here we show that all stable Brans-Dicke theories have 
this property. Setting Mpi = (SttG)"^/^ = 1, the Brans- 
Dicke action in Jordan frame takes the form 




{daf + C^[g]], (Bl) 



where Wbd is the Brans-Dicke parameter. In this frame 
the matter Lagrangian, C^n, is independent of a. 

For the study of perturbations, however, it is conve- 
nient to work instead in Einstein frame, obtained through 
the conformal transformation: 



5e 



6 ^'^QfJ.i^ : 



(B2) 



where in the last step we have introduced = 1/ (6 -|- 
4wbd)- Moreover, the field redefinition a = e^^^ is such 



that <j) is canonically normalized in Einstein frame: 

(B3) 

In the absence of anisotropic stress in the matter sec- 
tor, one can set = —"^e in Einstein frame. From the 
conformal transformation (B2) we can therefore read off 
in Jordan frame, to linear order. 



(B4) 



With the matter described by dust, it is straightfor- 
ward to show that the cosmological equations of motion 
allow for a scaling solution 



3 + 2/3^ 



4/3 



logtE, 



(B5) 



^ ' " 3 + 2/32 

where the subscripts make explicit that we are working 
in Einstein frame. Assuming adiabatic initial conditions, 
conservation of ( on super-horizon scales implies the fol- 
lowing growing-mode solutions 



S(l> = C— [ diEaE = C- 



2/3^-^3 ^ 
2/32 -h 5 ' 



2/32 + 5 ■ 



Substituting into (B4), we obtain 
4/32 



2/32 + 3 3a;BD + 5 



(B6) 



(B7) 



Absence of ghosts requires that 3 -|- 2wbd > 0, and thus 
9\k<^aH > for all stable Brans-Dicke theories, as we 
were hoping to show. 
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